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Motivated by a recent experiment with optical lattices that has realized a ladder geometry with 
an effective magnetic field (Atala et al., Nature Physics 10 , 588 (2014)), we study the dynamics of 
bosons on a tight-binding two-leg ladder with complex hopping amplitudes. This system displays 
a quantum phase transition even without interactions. We study the non-equilibrium dynamics 
without and with interactions, in the presence of a harmonic trapping potential. In particular we 
consider dynamics induced by quenches of the trapping potential and of the magnitude of the rung 
hopping. We present a striking ’’slowing down” effect in the collective mode dynamics near the phase 
transition. This manifestation of a slowing down phenomenon near a quantum phase transition can 
be visualized unusually directly: the collective mode dynamics can be followed experimentally in 
real time and real space by imaging the atomic cloud. 


I. INTRODUCTION 

The wavefunction of a charged quantum particle mov¬ 
ing in a magnetic held along a closed trajectory acquires 
a complex phase determined by the magnetic hux pen¬ 
etrating the enclosed area. An artihcial magnetic held 
can be designed for neutral particles by arranging for 
the wavefunction to acquire the appropriate amount of 
phase when the particle moves in space. This manner 
of creating artihcial gauge helds mm have been realized 
with photons [3] and in several cold atom setups, both in 
the continuum [11] and in optical lattices |7HI3|. In an 
optical lattice, a magnetic held can be mimicked using 
effective complex hopping amplitudes. This is done for 
example using laser-assisted tunneling. 

In Ref. [12], an effective magnetic held was created 
for bosons on a two-legged ladder lattice. In a tight- 
binding description, this can be described by adjusting 
the phases of the hopping amplitudes so that a hopping 
around one plaquette increases the phase by </>, which is 
then regarded as the hux per plaquette. In the Landau 
gauge this involves real hopping matrix elements along 
the legs and complex inter-leg coupling, with the phase 
of successive rung phases dihering by 4>. The experiment 
m explored the quantum phase transition between small 
and large magnitudes of the rung coupling, using hux 
4> = ^ per plaquette. The phase transition is between a 
vortex phase at small rung hoppings, where the ground 
state shows current patterns that can be described as 
a sequence of vortices along the ladder, and a so-called 
Meissner phase at large rung couplings. In the Meiss¬ 
ner phase, the magnetic held is “expelled” from the lad¬ 
der, so that the currents in the ground state run only 
around the boundary of the system, i.e, chiral currents 
along the two legs in opposite directions. The name de¬ 
rives from analogy to the physics of type-II superconduc¬ 
tors in a magnetic held: A type-II superconductor expels 
small helds (Meissner phase) but sustains a vortex lat¬ 
tice above a critical held. In addition to the characteristic 
current patterns, the vortex phase is also characterized 


by a modulated density prohle and a degeneracy of the 
single-particle spectrum. 

This phase transition was hrst studied theoretically in 
the context of Josephson ladders [HHii]- Several theo¬ 
retical studies of bosons in ladders subject to a magnetic 
held have appeared in the cold atom context, both be¬ 
fore and after the experimental realization of Ref. m- 
Ref. [T7], which is closely related to the experiment [T^ . 
studied the phase transition for non-interacting bosons. 
Several studies have appeared for interacting cases (Bose- 
Hubbard or hard-core bosons), using a variety of methods 
[THll25j . The predominant interest in these studies, ex¬ 
cepting [25], have been phase transitions and other equi¬ 
librium properties in uniform systems. 

In this work, we consider the non-equilibrium dynam¬ 
ics of bosons in the presence of an additional harmonic 
trapping potential. Cold atomic systems are particu¬ 
larly suitable for tracking non-equilibrium evolution in 
real time [MIIIT], and this possibility has generated an 
explosion of theoretical interest in the unitary dynamics 
of isolated quantum many-body systems [28] . Now that 
artihcial gauge helds have been realized in several setups, 
non-equilibrium evolution of such systems is a promising 
and exciting direction of future exploration. Here, we 
present a non-equilibrium study of one such system. We 
focus on trapped systems: The experiment m, like most 
cold-atom experiments, was performed in the presence of 
a trap. A many-body system in a trap has collective dy¬ 
namical modes which have no analog in uniform systems, 
the simplest being dipole and breathing modes. Because 
of the omnipresence of trapping potentials in cold-atom 
experiments, dipole and breathing mode dynamics are 
pervasive across the held: These modes have been inten¬ 
sively studied and used for diagnostic purposes since the 
beginning of this held, and continue to generate interest 
today [IZIdHin]. 

We study quenches (sudden changes) of the position 
and strength of the harmonic potential, which generates 
dipole and breathing mode dynamics respectively. For 
both non-interacting and interacting systems, we identify 
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in the spectrum the eigenstates associated with the dif¬ 
ferent modes, highlighting differences between the Meiss¬ 
ner and vortex phases and especially the behavior at the 
transition. The eigenstates associated with the modes 
are (very nearly) doubly degenerate in the vortex phase, 
but not in the Meissner phase. At the phase transition, 
the eigenstates rearrange their roles in order to accom¬ 
modate this difference. 

Close to the phase transition we observe a slowing 
down of collective mode oscillations — in the non¬ 
interacting case the dipole frequency even vanishes at 
the transition. This may be regarded as an analog of 
‘critical slowing down’ of dynamics near classical phase 
transitions. The effect found here should be experimen¬ 
tally detectable through direct imaging of the oscillat¬ 
ing size or position of the excited cloud. To the best of 
our knowledge, this is the only known visualization pos¬ 
sibility of a real-time slowing-down effect at a quantum 
phase transition. We attribute this effect to vanishing (or 
minima) of relevant energy gaps in the low energy spec¬ 
trum at the critical point. We also show that the effect 
persists in interacting systems, using exact diagonaliza- 
tion studies of the eigenspectrum and dynamics in multi¬ 
boson systems, and using the time-dependent Gutzwiller 
approximation |41) . This phenomenon provides an un¬ 
usual (real-space and real-time) manifestation of the re¬ 
lationship between vanishing energy scales and diverg¬ 
ing time scales near quantum phase transitions. The 
present results thus supplement, on the one hand, the 
various frequency-space consequences of this relationship 
[22], and on the other hand, the indirect manifestation 
of quantum-critical slowing down appearing in Kibble- 
Zurek-like scaling relationships for ramp (‘slow quench’) 
dynamics [2H]- 

We consider also quenches of the inter-leg hopping 
magnitude. This generates breathing modes, similar to, 
e.g,, interaction quenches in a trapped bosonic system 
without a magnetic field. We show that the difference 
of eigenstate arrangements in the two phases manifests 
itself in the particular combination of eigenstates which 
are excited, in a manner peculiar to this system. 

The paper is organized as follows. In Sec. jn] we in¬ 
troduce the Hamiltonian and notations. In Sec. inii we 
describe the low energy spectrum. We identify a relevant 
symmetry of the Hamiltonian and describe degeneracies 
of eigenstates. For context, we contrast with a pure one¬ 
dimensional system (e.g., single tight-binding chain) and 
connect excited states to collective modes. In Sec. IIVI we 
study quenches of the trapping potential. This section 
contains the most striking result of the paper, namely 
critical slowing down of collective modes close to the 
quantum critical point. We present and analyze exact 
time evolution results for systems up to four bosons as 
well as mean-field (Gutzwiller) results for seven particles. 
We connect the observed dynamical slowing down to the 
energy spectrum. In Sec. jV] we show results for inter¬ 
leg hopping magnitude quenches. Sec. |VI| summarizes 
the results and presents open problems. In Appendix fAj 


we present the effects of the geometry of finite ladders 
(boundary conditions, presence of trap) on the lower part 
of the energy spectrum. 


II. SYSTEM AND HAMILTONIAN 

The Hamiltonian we consider is 

H = -jy, e 

I (7 — L^R 

I 

Z,cr Z,cr 

Here bi-u, b\.^ are bosonic operators for the site on rung 
I {I = 1... L) and leg a = L,R, ni-a = bj.^b^a, 4> is the 
magnetic flux per plaquette, J is the intra-leg and 
is the inter-leg hopping matrix element, Iq is the center 
of the trap, and fctr is the trapping strength. We take 
the trap to be centered at one of the rungs, i.e., Iq is an 
integer. 

It is possible to use a gauge in which the complex phase 
is on the leg hoppings instead of the rung hoppings, as 
in Refs. [Mini El [SUM! • If the flux (j) per plaquette is 
the same, physical properties are unchanged. 

When the density is low enough, we can approximate 
the cosine dispersion of a lattice particle by a quadratic 
dispersion (in the leg direction). This “effective mass ap¬ 
proximation” ascribes a continuum mass to lattice parti¬ 
cles, m* = ^, so that we can relate our trapping strength 
fctr to the trapping frequency of a continuum trapping po¬ 
tential ujq = \/2Jktr- We use this expression 

as the definition of the trapping frequency uiq. 

Fig. s a) shows a sketch of the magnetic ladder sys¬ 
tem. The phase transition between Vortex and Meiss¬ 
ner phases can be tuned either by magnetic flux (j) or 
by inter-leg hopping K m- Here we fix </) = ^ as in 
the experimental realization and change K. The value of 
K at the critical point is Kc = for ^ = |. We use 
h = 1, J = 1, measure energy in units of J and use lattice 
spacing a = 1. 

In the absence of a trap or interactions (/ctr = U = 0), 
the phase transition can be understood by considering 
the single-particle dispersion as a function of momentum 
[H EH ESI, which is well-defined for translationally in¬ 
variant systems, e.g., for infinite systems and for periodic 
boundary conditions (pbc) with L = (integer n). 

For K < ATc, the dispersion has two minima, at nonzero 
momenta of opposite signs. Thus in the vortex phase 
the ground state is doubly degenerate and has finite mo¬ 
mentum. When translational invariance is broken, such 
as with open boundary conditions (obc) or with a trap, 
momentum is no longer a good quantum number, but an 
approximate degeneracy persists, and the density profile 
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shows oscillations with a wavelength corresponding to the 
inverse of the above-mentioned finite-momentum value. 
In the Meissner phase {K > K^), the dispersion has a sin¬ 
gle minimum at zero momentum; the eigenstates are not 
degenerate, and the density profile in obc and trapped 
cases has no modulations. 


III. LOWEST EIGENSTATES OF TRAPPED 
SYSTEM 


In this section we describe the low-energy spectrum 
of the magnetic ladder with an additional harmonic po¬ 
tential. For harmonically trapped systems, the low¬ 
est excited states have simple dynamical interpretations: 
These states correspond to various collective excitations. 
Before treating the magnetic ladder, we first review the 
case of a pure ID system (either continuum or a sin¬ 
gle tight-binding chain with small fillings) in a harmonic 
trap. 


A. Review: one-dimensional system 


The ID bosonic Hamiltonian with harmonic trapping 
potential is reflection-symmetric with respect to the trap 
center. The eigenstates are alternately symmetric and 
antisymmetric under reflection. In the single-particle 
case, these eigenstates are the familiar harmonic oscil¬ 
lator eigenstates, and are characterized by an increasing 
number of nodes, e.g., the ground state is nodeless, and 
there are one and two nodes in the first and second ex¬ 
cited eigenfunctions respectively. 

The lowest two excited states are associated respec¬ 
tively to dipole mode and breathing mode excitations. 
The dipole mode involves oscillations of the center-of- 
mass position. It can be excited by shifting the position 
of the trap minimum, i.e., a trap position quench. The 
breathing mode involves oscillations of the width of the 
density distribution, i.e., expansion and contraction of 
the cloud. It can be excited by a quench of either the trap 
strength, or (in interacting systems) through a quench of 
the interaction strength. 

The frequency of the collective mode oscillations is 
equal to the excitation energy of the corresponding eigen¬ 
states above the ground state, i.e., ujum = Edm — Eq 
and ujbm = Ebm — Eq for the dipole and breathing 
mode respectively. The dipole mode frequency ujbm is 
independent of interaction strength, and is locked to the 
trapping frequency (Kohn’s theorem |4 ()] ). In a non¬ 
interacting system, the spectrum is equally spaced, so 
that the breathing mode frequency ujbm is twice the trap¬ 
ping frequency, i.e., the ratio between the lowest collec¬ 
tive modes is ujbm/^dm = 2. For an interacting system, 
this ratio depends on U and has a minimum at inter¬ 
mediate U, which for large systems is ujbm /ujdm = 

[SI1IS1I37]. 



FIG. 1. (a) Sketch of the magnetic tight-binding ladder 

with complex inter-leg hoppings. Hopping around a plaque- 
tte increases the phase by (f>. (b) Low energy excitation 

spectrum of the trapped magnetic ladder with ^ = f for 
N = 1, L = 500, fctr = 10““^. The eigenstates are nearly 
degenerate for K < Kc (vortex phase) and non-degenerate 
for K > Kc (Meissner phase), (c) Density distributions of 
four lowest eigenstates for N = 1,L = 120, fctr = 10“® and 
K = {0.5;!} < Kc and K = {1.5; 2} > Kc on a leg. For 
K < Kc the degenerate eigenstate pairs have identical den¬ 
sity distributions. The first degenerate pair have gaussian-like 
shape and the second pair corresponds to the dipole mode. 
For K > Kc the states are rearranged: the second eigenstate 
takes the shape of a dipole mode state and the third eigen¬ 
state becomes a breathing mode state. The labels mark the 
corresponding eigenenergy in panel (b), and the ©, 0 signs 
indicate parity under pseudo-inversion. 


B. Pseudo-inversion symmetry of the magnetic 
ladder Hamiltonian 

Because of the space-dependent complex phase in the 
rung hoppings, the magnetic ladder Hamiltonian 0 is 
not symmetric under reflection about the trap center. 
However, if the trap is centered on the central rung of 
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the ladder {Iq = ) or if the ladder is infinite, then the 

Hamiltonian does possess a pseudo-inversion symmetry. 
The symmetry transformation involves reflection about 
the trap center and exchange of legs (reflection around 
rung center), together with a space-independent phase 
factor. Formally the symmetry transformation leaving 
the Hamiltonian invariant is Each 

eigenstate of the Hamiltonian is either odd or even under 
this pseudo-inversion transformation. 

Because simple reflection around the central rung is 
not a symmetry of the Hamiltonian, the eigenstates are 
not required to have definite parity under such reflec¬ 
tion. There is no straightforward correspondence be¬ 
tween the density profile (Gaussian-like, dipole-shaped 
etc) and parity, because the well-defined parity in this 
system corresponds to pseudo-inversion and not reflec¬ 
tion. 

Note that, in order to use this symmetry, we do not 
need the trap center to be at the ladder center as long as 
we are considering eigenstates that are localized enough 
that they do not feel the edges of the chain: such eigen¬ 
states will possess well-defined parity with respect to in¬ 
version around the center of the trap. The symmetry is 
also present if the trap center is between two rungs, i.e., 
lo is half-integer. 


C. Single-particle eigenspectrum of trapped 
magnetic ladder 

We now describe the eigenspectrum of the magnetic 
ladder with a trap for a single particle, for varying rung 
hopping K. The non-interacting multi-particle spectrum 
can be constructed by filling up the single-particle states. 
The interacting multi-particle case will be treated briefly 
in a later section. 

Figure [^b) shows the lowest few states of the single¬ 
particle spectrum, for a weak trap. Clearly, the tran¬ 
sition (between degenerate states at small K and non¬ 
degenerate states at large K) survives almost unaltered 
in the presence of a trap: the low-lying states are (nearly) 
degenerate for small K and non-degenerate for large K. 
The point K,. where the two lowest states split is a bit 
smaller than in the uniform case, but —>■ y/2 for de¬ 

creasing trapping strength, fctr 0. The excitation en¬ 
ergy of excited states have a minima near K^- 

For weak traps (fctr ^ 10“"^), the eigenenergy pairs 
for K < Kc are extremely close, the energy differences 
between the lowest two states being well below double¬ 
precision accuracy (details in Appendix 0. Thus for 
most purposes they can be regarded as degenerate pairs, 
and we will refer to the pairs as such, even though the 
degeneracy is not strictly speaking exact. Each such ‘de¬ 
generate’ pair includes one even and one odd eigenstate 
under the pseudo-inversion symmetry introduced above. 

At AT —>• 0 the two legs of the ladder are two separated 
tight-binding chains. Energy eigenstates are thus truly 
doubly degenerate. The lowest eigenstates are wo-spaced. 


in contrast to the uniform case. At AT —>■ oo the magnetic 
ladder behaves as a single tight-binding chain with each 
effective ‘site’ composed out of a rung, resulting in an 
effective hopping amplitude J* = cos(|)J = J/-\/2 ^0] 
and effective trapping frequency Wq = 2“Ja;o. In this 
limit the states are non-degenerate and Wg-spaced. 

Figure [^c) shows plots of the density distributions 
of the four lowest eigenstates at different AT along one 
leg. The density distribution is equal on both legs. For 
AT > ATc (lower two panels) the low energy spectrum 
contains non-degenerate states whose density profiles are 
familiar from single-chain eigenstates: they have an in¬ 
creasing number of nodes for increasing energy, and cor¬ 
respond to a single gaussian eigenstate, a dipole eigen¬ 
state, a breathing eigenstate, etc. For AT < ATc, the 
low-lying eigenstates all have density modulations, as ex¬ 
pected for the vortex phase [12]. The first two states 
(belonging to the ground state manifold) have Gaussian- 
like density profiles (with additional vortex modulation). 
The second pair have one ‘node’ each in the density dis¬ 
tribution, i.e., both have the form of a dipole eigenstate. 
Thus the gaussian, dipole, breathing,... eigenstates each 
come in degenerate pairs. Each pair has one inversion- 
symmetric and an inversion-antisymmetric eigenstate. 

At the transition point, the states change their na¬ 
ture. The ground state manifold splits into a gaussian- 
shaped ground state and a dipole-mode-shaped excited 
state. The other degenerate manifolds similarly change 
shape: the lower of the second degenerate pair changes 
from a dipole profile to a breathing mode profile having 
two ‘nodes’. 

The lifting of the degeneracy of the ground state at the 
critical point ATc, and the lowering of the other excita¬ 
tion energies around K = K^, have intriguing dynamical 
consequences, detailed in the next section. The dipole 
and breathing frequency decrease dramatically while ap¬ 
proaching Kc- 


IV. TRAP QUENCHES: CRITICAL SLOWING 
DOWN 

In this section we present ‘critical slowing down’ dy¬ 
namics of the collective modes, both for single- and multi¬ 
boson systems (for N > 1). We excite collective modes 
by quenching the trapping potential, i.e. trapping po¬ 
sition and trapping strength for dipole and breathing 
modes respectively. We calculate time dynamics of first 
and second moments of the density distribution, i.e., 
center-of-mass position and ‘width’ respectively. Close 
to the Meissner-vortex phase transition in a trapped sys¬ 
tem the frequencies of collective modes are significantly 
suppressed, i.e., the system slows down. 
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FIG. 2. (Color online) Dipole mode oscillation after a quench (sudden shift) of trap position, (a) Real time dynamics of 
center-of-mass in both legs (L,R), after a trap shift of do = 1 for a system with L = 300, fctr = 0.001 for several K values. The 
slowing down effect is clearly visible: Oscillation period is largest close to the critical point, (al) Exact result for = 1. (a2) 
The time-evolution is calculated using Gutzwiller mean-field approximation for = 7,17 = 1. (b) Increase of dipole mode 
oscillation period close to the critical point at Kc- (bl) Exact results: non-interacting {N = 1) and interacting {N = 4) systems 
compared. At small K and large K, the frequencies become identical in trap units (Kohn theorem). (b2) Same effect is seen 
with mean-field dynamics. 


A. Trap position quenches: dipole mode dynamics 

In the Figure l^al) we plot oscillations of the center- 
of-mass position after excitation of the dipole mode for 
a single particle {N = 1) at different inter-leg hopping 
magnitudes K. The dipole mode is excited by initializing 
the system in the inversion-symmetric ground state and 
then shifting the center of the potential by one site, do = 
1 . 

The frequency of the dipole mode oscillation is strongly 
dependent on K. The dipole frequency is wg at A = 0 
and Wg = 2“Jwg at A —)■ oo. In between the fre¬ 
quency has a pronounced minimum around Kc, as vis¬ 
ible through the several times larger periods at K « 1.2. 
The dipole frequency follows closely the excitation en¬ 
ergy of the first excited state for K^Kc and the ex¬ 
citation energy of the first excited degenerate pair for 
K<^Kc. This suggests that the frequency could vanish 
a.t K = Kc, since the lowest two states merge as is ap¬ 
proached from above. However, in the transition region, 
the dynamics becomes multi-frequency, because multiple 
eigenstates get excited. (This is discussed below in the 
multi-particle context. Figure]^) 

The same slowing-down effect is observed in multi¬ 
boson systems, also in the presence of interactions. In 
the panels in Fig. [^a2), we plot mean-field results for 
N = 7 bosons, U = 1, obtained using the time-dependent 
Gutzwiller approximation |41] . The Gutzwiller approxi¬ 
mation is of course not exact, but this calculation demon¬ 


strates that the slowing down effect is generic and ap¬ 
pears also with non-negligible interactions. 

In Figure [^bl,b2) we show the dipole oscillation pe¬ 
riod for A = 1, for A = 4 interacting bosons obtained 
using exact numerical diagonalization, and for N = 7 
interacting bosons using the Gutzwiller approximation. 

The dynamics shown in Fig. is after a trap shift by a 
single lattice site, dg = 1. For stronger quenches (dg > 1) 
the time evolution is overall similar and the same slowing 
down effect is seen. For larger dg, more states are excited 
resulting in stronger beating effects. 

Dipole excitations of harmonically trapped one¬ 
dimensional system are subject to Kohn’s theorem |40j . 
In the present case, only the K = 0 and K = oo lim¬ 
its can be regarded as effectively one-dimensional sys¬ 
tems. In accordance we see in Fig. [^bl) that for small 
and large values of K, the oscillation period is equal for 
non-interacting and interacting system at 17 = 1, when 
measured in units of the trapping frequency. There is no 
analogy of Kohn’s theorem for a finite K. 


B. Trap position quenches: excitation of 
eigenstates 

The dynamics after a quench is determined by the 
eigenstates of the final Hamiltonian which are excited 
through the quench, i.e., the distribution of overlaps 
of the eigenstates of post-quench Hamiltonian with the 
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FIG. 3. (Color online) Excitation energy spectrum of the 
trapped magnetic ladder with N = 2 bosons, L = 50, k = 
0.001 and (a) t/ = 0 and (b) U = 1. The line thicknesses 
logarithmically encode overlaps \{4>o\'tpi)\'^ of eigenstates with 
initial ground state of the Hamiltonian with a trap shift by 
do = 0.5, i.e. after dipole mode excitation. The thick, colored 
lines correspond to states with an overlap larger than 10~^. 
Dashed lines correspond to states with overlap smaller than 
10 ~^. 


initial state, \{(j)o\'ijji)\'^. The energy difference between 
states having substantial overlap are the dominant fre¬ 
quencies in the time evolution after a quench. 

Fig. [3] shows the energy spectrum for N = 2 bosons 
without and with interactions. The overlap distribution 
after a trap position quench at each value of K is encoded 
in the thickness of lines. 

Without interactions {U = 0, Fig. ia)) , the > 1 
spectrum has the same qualitative features as for TV = 1, 
except for additional degeneracies at K < Kc and at 
large K^K^, and visibly larger number of states around 
K Ki Kc- The larger number of states and greater degen¬ 
eracies can be understood as the number of possibilities 
for multiple bosons to occupy the single-particle states. 
Away from the critical region, the effect of a trap position 
quench is to excite mainly the new ground state and the 
first excited state (dipole eigenstate), so that the lowest 
energy gap sets the dipole frequency. In the AT « ATc 
region, several states are seen to be excited, leading to 
multi-frequency dynamics reported already in Fig. [^a). 

Interactions lift degeneracies in the energy spectrum, 
Fig.§b). For K < Kc, this leads to the situation, per¬ 
haps surprising at first sight, that the dipole-mode eigen¬ 
state is not the lowest excited state. This does not occur 
in pure one-dimensional (single-chain) systems [31] [33] 
and is a novel feature of ladder geometries. The exci¬ 
tation energies of the lowest excited states have a min¬ 
imum at A' « K(. but there is no proper gap closing 
as one approaches from the AT > A'c side, in contrast 
to the U = 0 case. As in the U = 0 case, the over¬ 
lap distribution after a trap quench is dominated by the 
ground state and dipole eigenstate away from the critical 
region, while there are multiple eigenstates excited in the 
AT « Kc region. For K < Kc, the excited states which 
appear below the dipole-mode state are not excited. The 
dipole-mode eigenenergy difference with the ground state 
energy is wq at AT = 0 and Wq = 2“4Wo at AT —>■ oo, just 
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FIG. 4. (Golor online) Breathing mode oscillation after a 
trap quench for various K, shown using the second moment 
of the density distribution for L = 300, ktr = 0.00105 —>■ 
0.001. (a) Exact result for N = 1. (b) Gutzwiller mean-field 
approximation for N = 7,U = 1. (c) Ratio of the breathing 
and dipole mode frequencies for N = 1, L — 300, kn = 
0.0001. The two excitations are obtained using a do = 1 
quench (dipole mode) and a fetr = 0.00011 —0.0001 quench 
(breathing mode). 


as in the U = 0 case, in accordance with Kohn’s theo¬ 
rem according to which this energy difference should be 
[/-independent. 


C. Trap strength quenches: breathing mode 

A similar slowing-down effect is visible for the breath¬ 
ing mode oscillations, in Fig. |^a,b). The frequency re¬ 
duction is smaller compared to the dipole mode case. 

In Fig. I^c) we plot the ratio of the energy gaps corre¬ 
sponding to breathing and dipole mode oscillations, cal¬ 
culated from the energy spectrum as a function of inter¬ 
leg hopping K. The energy gap used is that between the 
ground state and the next most strongly excited state in 
a trap strength quench (for the breathing mode) and in a 
trap position quench (for the dipole mode). Away from 
the critical point, the ratio is equal to two, as expected for 
a single particle. Close to the critical point, the behav¬ 
ior is completely different from non-magnetic system; the 
divergence-like feature signals that the dipole frequency 
nearly vanishes in this region. 


V. INTER-LEG HOPPING QUENCHES 

In this section we treat quenches of the inter-leg hop¬ 
ping magnitude K. In Fig. [^ we fix the initial Ki to be 
in the vortex phase, Ki < and show overlaps with 
eigenstates of the final Hamiltonian for various Kf. We 
choose the initial state to be one of the two states of the 
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FIG. 5. (Color online) Excitations after a IF-quench. We 
show overlaps of the eigenstates with an initial state, encoded 
in the thickness of the solid lines. The initial state is one state 
of the ground state manifold with Ki = 0.5, N = 1, L = 200, 
fctr = 0.001 with (a) and (b) odd parity under inversion. Only 
low energy states with the same symmetry as the initial state 
are excited. 


ground state manifold for Ki < K^-, i.e. either symmetric 
or anti-symmetric with respect to the pseudo-inversion 
transformation introduced in subsection IIIIB1 Note that 
both these states have gaussian-like density profiles with 
additional density modulations (Fig.j^c), states 1,11,1,2). 

For Kf < Kc, the quench chooses from each degen¬ 
erate pair the state with the same symmetry as the ini¬ 
tial state, and one of each degenerate pair has nonzero 
overlap. The highest overlap is with the ground state, 
and the overlaps get progressively smaller with increasing 
eigenenergy. For Kj > Kc, alternate eigenstates have the 
same symmetry, so that the nonzero overlaps are either 
only with odd states, or only with even states. The two 
cases are shown in the two panels of Figure The situa¬ 
tion on the right panel may appear counter-intuitive, as 
the eigenstates excited have density profiles that appear 
odd, which for a pure one-dimensional system one would 
not expect starting from a state with a gaussian profile. 
However, the more complicated symmetry of the present 
case allows odd-symmetry states to have gaussian den¬ 
sity profile and even-symmetry states to have dipole-like 
density profile. 

The temporal dynamics after a K quench involves 
complicated multi-frequency oscillations including both 
breathing and dipole-like components in the densities of 
both legs (not shown). The total center-of-mass stays 
fixed but there are out-of-phase dipole oscillations in the 
two legs. 


VI. SUMMARY AND DISCUSSION 

In this work, we have provided a theoretical study of 
some non-equilibrium dynamics of bosons in the mag¬ 
netic ladder system recently realized in experiment |12j . 
focusing in particular on effects due to the presence of a 
harmonic trapping potential, which is almost universally 
present in such experiments. Previous theory studies 


have focused heavily on the uniform system and equilib¬ 
rium properties. Other than this work, non-equilibrium 
non-uniform dynamics was also studied in Ref. [25) . 
where a nonlinear version of the system was studied in the 
presence of disorder. This may be regarded as a mean- 
field approximation to a many-boson system with on-site 
interactions. 

By studying trap quenches, we have presented the dy¬ 
namics and spectral structures related to dipole modes 
and breathing modes. The most striking result is a 
slowing-down effect in these modes near the critical point, 
an effect that should be readily observable in the exper¬ 
imental setup. We are unaware of any other example of 
a possibility to visualize quantum critical dynamics in 
real-space oscillation phenomena. Critical slowing down 
in real-time dynamics is of course well-known for clas¬ 
sical phase transitions [43], where the phenomenon can 
be dramatically visualized, e.g., through the response of 
a liquid sample near a liquid-gas transition point after 
being shaken. 

We have stressed in particular the connection of spe¬ 
cific dynamical phenomena to structures in the eigen- 
spectrum. This correspondence, and the fact that note¬ 
worthy spectral structures can rise to unusual dynam¬ 
ical phenomena, is increasingly appreciated and explic¬ 
itly studied in the recent non-equilibrium literature, e.g., 
in the context of self-trapping and related phenomena 
in Bose-Hubbard dimers and ladders m 113, in the 
study of (repulsive) binding of lattice particles |46| , col¬ 
lective modes of trapped systems imiMj, dynamics in 
spin chains [17], and various other phenomena [451150j . 
Our results connecting critical slowing down of collective 
modes to minima in the excitation gaps (or splitting of 
near-degenerate eigenstate pairs) can be seen as an ad¬ 
dition to this theme. In addition, in this work we have 
examined the eigenstates of the final Hamiltonian that 
are excited in a quench, i.e., the distribution of overlaps 
of the final eigenstates with the initial state as a function 
of energy. This overlap distribution excited in a quench, 
and its continuous version, the work distribution, has 
been the focus of much recent attention [5TH55] . 

A peculiar feature of the present system is that 
gaussian-like and dipole-like shapes do not necessary cor¬ 
respond to even and odd parity respectively, in contrast 
to the experience obtained from studying purely one¬ 
dimensional trapped systems or the harmonic oscillator 
eigenstates. The reason is that the conserved parity in 
this case does not relate to a simple reflection around 
the center rung, but rather a pseudo-inversion involv¬ 
ing reflections along both leg and rung directions (plus a 
possible phase shift). This allows somewhat unexpected 
situations like the first two (near-degenerate) eigenstates 
both having gaussian-like profile, the second pair hav¬ 
ing dipole-like profile, etc., in the vortex phase. At the 
transition, the states rearrange and change nature, and 
in the Meissner phase we get a more familiar sequence of 
a single ground state with gaussian-like density profile, 
only the first excited eigenstate with dipole-like profile. 













FIG. 6. (Color online) Energy gap between lowest two single¬ 
particle eigenstates, obtained using double-precision numer¬ 
ics. (a) Open boundary conditions, various sizes. The gap in 
the vortex phase is seen to be numerically significant, 

(b) Trapped system with various ku- Data is obtained with 
L = 4000, for which the trapped cloud does not reach the 
edges up to numerical accuracy. The gap is seen to reach 
“numerical zero”, ~ 10“^®. 


der systems which focus on uniform systems. To provide 
context, we describe here some effects of the geometry of 
the system, by contrasting spectral features of periodic, 
open-boundary and trapped systems. 

Periodic boundary conditions — When L is a mul¬ 
tiple of 27r/0 = 4, a system with pbc is uniform, and 
momentum along the leg direction is a good quantum 
number. In such a case diagonalization can be done ana¬ 
lytically in the absence of interactions, by Bogoliubov ro¬ 
tation in momentum space. The resulting single-particle 
spectrum has two branches, which we find to be 

£1(1) = -2JAq - 2KCq , £1(2) = -2JBq + 2KCq (Al) 

with 

Aq = cos q sin^ 0 -I- cos{q + 4>) cos^ 0, 

Bq = cos ((7 -I- (j)) sin^ 0 -I- cos q cos^ 0, 


and so on. 

We have also shown how interactions modify the spec¬ 
trum for a multi-boson trapped system. Interactions lead 
to the peculiar effect that, in the Meissner phase, the 
dipole eigenstates are no longer the first excited state or 
first excited pair; instead there are eigenstates between 
the ground state and the dipole eigenstate. This is a sit¬ 
uation that does not occur in pure ID chains mmm- 

The present work raises a number of open questions. 
The very nearly perfect degeneracy in trapped finite sys¬ 
tems remains unexplained. As we show in the Appendix 
(Fig. §, open boundary conditions cause degeneracy¬ 
breaking gaps in the Meissner phase that are many orders 
of magnitude larger than in even more tightly confined 
harmonic-trapped cases. Analytic solutions for the obc 
and trapped cases would be helpful in this regard but 
are currently not known. Physically, one can speculate 
that harmonic trapping is a less violent manner of break¬ 
ing translation symmetry than open ends of the ladder. 
More generally, this is one of the first studies of non¬ 
equilibrium dynamics in the magnetic ladder, treating 
the most obvious quench protocols and the correspond¬ 
ing parts of the eigenspectrum. We expect many other 
interesting dynamical phenomena, e.g,. corresponding 
to aspects of the eigenspectrum that have not yet been 
explored, or involving interesting types of current circu¬ 
lation dynamics. 
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Appendix A: Geometry 

In this work we concentrate on trapped systems, in 
contrast to most of the literature on the magnetic lad¬ 


Cq = sin 0 cos 0, 

and 

^ 1 K \ 

0 = - arctan —-^--r-. 

2 \ J(cos(g -I- - cos(g))y 

Momenta take the values q = with an integer n = 
-£/2,-£/2 + l,...,£/2- 1. 

When L is an even multiple of 4, i.e., a multiple of 
8, momenta values can be shifted by cj)/2: k = q — (/)/2 
recovering the expression appearing previously in the lit¬ 
erature [ninnii^: 


E± = —2J cos k cos ^ ± Y 4 sin^ k sin^ ^ -I- . (A2) 

The spectrum contains two minima for K < (hence 

the exact degeneracy) and one minimum for K > 
(hence the single ground state). Note that this is correct 
only for L = 8n with integer n. For £ = 4n with an odd 
n energy spectrum contains doubly degenerate states fro 
all K. 

When L is not a multiple of 4, we are unaware of an an¬ 
alytical solution for pbc without any twist in the bound¬ 
ary condition. Of course, expression ( |A2[ ) was previously 
used with large sizes in mind, where the boundary condi¬ 
tions become insignificant. Our purpose here is to present 
explicitly the effects of geometry. 

When momentum is well-defined (£ = 4n), almost all 
states are degenerate. The only exceptions are the k = 0 
states for even n {L = 8n). At small K, the full spec¬ 
trum (not shown) shows many true level crossings be¬ 
tween such degenerate pairs. 

Open boundary conditions — With obc, momentum 
is not a good quantum number. Level crossings are there¬ 
fore avoided, so that the full spectrum looks visually quite 
different from the pbc case. Of course, for large enough 
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system sizes, the low-energy part of the spectrum looks 
similar to the pbc case. For K < Kc, the energy levels 
are grouped in pairs, each such pair undergo a succes¬ 
sion of mutual avoided level crossings as a function of K. 
With increasing L, the gap within each pair decreases, as 
shown in Fig.j^a). This is consistent with the obc case in 
the large-L limit becoming similar to the pbc case (where 
the K < Kc pairs are degenerate). 

System with harmonic trap — The energy spectrum 
of the trapped system is similar to a system with obc; 


momentum is not a good quantum number. However, 
the trapped system displays extremely small gaps be¬ 
tween the degenerate pairs for K < Kc- In Fig. [^b) we 
show the lowest gap as obtained from double-precision 
computations. The pairs are ‘degenerate’ up to standard 
numerical accuracy already for fctr ~ 10“^. Numerical 
diagonalization with higher precision shows that the de¬ 
generacy is not exact and there is a minute but hnite gap 
for finite ktv At present a detailed explanation of this 
ultra-small gap is lacking, since we do not have analytical 
solutions for the trapped or obc cases. 
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